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ABSTRACT : Interferon regulatory factor 6 (IRF6) gene is a member of the IRF-family, and plays functionally diverse roles in the 
regulation of the immune system. In this report, the 13,720 bp porcine IRF6 genomic DNA structure was firstly identified with a 
putative IRF6 protein of 467 amino acids. Alignment and phylogenetic analysis of the porcine IRF6 amino acid sequences with their 
homologies to other species showed high identity (over 96%). Tissues expression of IRF6 mRNA was observed by RT-PCR, the results 
revealed IRF6 expressed widely in eight tissues. One SNP (HQ026023:1383 G>C) in exon7 and two SNPs (HQ026023:130 G>A; 232 
C>T) in the 5 promoter region of porcine IRF6 gene were demonstrated b y DNA sequencing analysis. A further analysis of SNP 
genotypes associated with immune traits including IFN-y and ILIO concentrations in serum was carried out in three pig populations 
including Large White, Landraces and Songliao Black pig (a Chinese indigenous breed). The results showed that the SNP 
(HQ026023:1383 G>C) was significantly associated with the level of IFN-y (d 20) in scruin (p = 0.038) and the ratio of IFN-y to ILIO (d 
20) in serum (p = 0.041); The other two SNPs (HQ026023:130 G>A; 232 C>T) were highly significantly associated with ILIO level in 
serum both at the day 20 (p = 0.005; p = 0.001) and the day 35 (p = 0.004; p = 0.006). Identification of the porcine IRF6 gene will help 
our further understanding of the molecular basis of the IFN regulation pathway in the porcine immune response. All these results should 
indicate that the IRF6 gene can be regarded as a molecular marker associated with the ILIO level in serum and used for genetic selection 
in the pig breeding. (Key Words : Pig, IRF6, Expression, Polymorphisms, Association Analysis) 



INTRODUCTION 

The interferon regulatory factor (IRF) family of 
transcription factors was initially found to be involved in 
the induction of genes that encode type I IFN and as 
essential regulators of the activation of immune cells with 
the discovery of pattern recognition receptors (PRRs), 
which function as a 'platform' that links innate and adaptive 
immune response (Akira et al., 2006). The IRF gene family 
consists of nine members encoding transcription factors that 
share a highly conserved helix-turn-helix DNA-binding 
domain and a less conserved protein-binding domain. IRFs 
play important roles in immune response to viral infection, 
cytokine signaling, cell growth regulation and 
hematopoietic development (Nguyen et al., 1997; Honda et 
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al., 2006) and most IRFs regulate the expression of IFN 
after viral infection. 1RF6 gene is one of nine members of 
IRF-family and has been identified in human, mouse, and 
zebrafish (Kondo et al., 2002; Ben et al., 2005). Mutations 
in IRF6 cause Van der Woud syndromes (VWS) and 
Popliteal pterygium syndrome (PPS) has been reported for 
human (Lees et al., 2000; Kondo et al., 2002). However, the 
function of IRF6 remains unknown in pigs. In order to gain 
more knowledge of the porcine IRF6 gene regulating IFN 
involved in resistance to viral infection, we firstly identified 
the genomic structure of porcine IRF6 gene, and 
characterized its molecular organization, detected mRNA 
expression patterns and polymorphisms. 

Cytokines are important mediators of the immune 
responses. The level of cytokines in serum such as IFN-y 
and ILIO, vary with the organic state of health. IFN-y is 
critical for iimate and adaptive inmiunity against viral and 
intracellular bacterial infections (Schoenbom and Wilson, 
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2007). ILIO has pleiotropic effects on immunoregulation 
and inflammation, it down-regulates the expression of Thl 
cytokines and MHC class II antigens (Redpath et al., 1999). 
Differences in levels of these immune traits and their ratios 
in serum among individuals under the same condition 
should provide evidence of genetic control. Considering the 
influences on IFN immune response and cytokines 
signaling pathway of IRFs genes, IRF6 could be an 
important candidate gene for porcine immune response and 
disease susceptibility, we also further performed an 
association analysis between the SNPs of the IRF6 gene and 
immune traits (IFN-y, IL-10 and IFN-y/IL-lO) in serum to 
estimate their possible effects in pigs. 

MATERIALS AND METHODS 

Animal and tissue samples 

There are 3 pig breeds in our study, including Landrace, 
Large White and Songliao Black (a Chinese indigenous 
breed). They were sampled from the experimental farm of 
the Institute of Animal Sciences, Chinese Academy of 
Agricultural Sciences, Beijing, China. All pigs were 
vaccinated with (Classical Swine Fever, CSF) live vaccine 
at the age of 21 d. The first blood samples were collected 
from each piglet one day before the vaccination (d 20), and 
2 wk after the vaccination, the second blood samples were 
collected (d 35). Ear tissue samples of all pigs were also 
collected for DNA extraction. Nine different tissues, 
including heart, liver, spleen, lung, kidney, skeletal muscle, 
lymph node, brain and thymus, were collected from three 
35-d-old Landrace pigs, then immediately frozen in liquid 
nitrogen and stored at -80°C for spatial expression analysis. 

Measurement of serum IFN-y and IL-10 concentration 

Serum samples from three pig populations were stored 
at -80°C. A commercial ELISA kit (Biosource, Carlsbad, 
California) was used to measure the immune traits 
including IFN-y and IL-10 concentration in each serum 
sample based on the standard instructions from 
manufacturer. The kit contains a antibody coated 96-well 
test plate, standards of known IFN-y/ILlO concentrations, 
standard diluent buffer, wash buffers, and a biotin 
monoclonal antibody specific to IFN-y/ILlO, a streptavidin- 
HRP diluent solution and a stop solution. All samples were 
randomly placed in each plate and were used to calculate 
IFN-y and IL-10 concentrations based on a standard curve. 

Total RNA extraction and genomic DNA isolation 

The total RNA was extracted from different tissues with 
a TRIzol reagent (Invitrogen, USA). In order to prevent 
serious contamination of genomic DNA, DNasel (Beijing 
Tiangen Biotechnology, China) treatment on the total RNA 



was carried out before first-strand cDNA synthesis. 
Genomic DNA was isolated from the ear tissue samples of 
Landrace, Large White and Songliao Black pigs using 
phenol/chloroform extraction and ethanol precipitation 
(Sambrook et al., 1989). 

RT-PCR analysis of IRF6 mRNA expression 

Based on the cloned mRNA sequence of IRF6 gene 
(GenBank Accession Number: AF327368), the reverse 
transcriptase-polymerase chain reaction amplification 
method was used to detect mRNA expression of the IRF6 
gene. The first strand cDNA was synthesized in the 
presence of 2 ^g total RNA, 0.5 oligo(dT)18, 200 nM 
dNTPs, 10 U RNAsin (Promega, Madison, WI, USA), 
IxM-MLV RT buffer, and 300 U M-MLV reverse 
transcriptase (Promega, USA) in a volume of 50 (xl at 40°C 
for 1 h. A specific primer setting for IRF6 gene mRNA 
expression analysis was designed to eliminate potential 
confounding results from genomic DNA contamination 
(Table 1). Another specific primer was also designed to 
amplify porcine house-keeping gene glyceraldehyde3- 
phosphate dehydrogenase (GAPDH) mRNA (GenBank 
accession number: AF017079) as a positive control 
(GAPDH-Fl 5 -TGAGACACGATGGTGAAGGT-3' and 
GAPDH-Rl 5'-GGCATTGCTGATGATCTTGA-3'). 

Genomic organization and SNP identification of the 
porcine IRF6 gene 

To characterize the genomic structure of the porcine 
IRF6 gene, the mRNA sequences of porcine IRF6 gene 
(GenBank Accession Number: AF327368) was used to 
search for porcine homologues with high throughput 
genomic sequences (HTGS) in the porcine genomic 
databases by BLAST (http://blast.ncbi.nlm.nih.gov) of 
NCBI GenBank database (http://www.ncbi.iilm.mh.gov/). 
One genomic DNA fragment (GenBank accession Number: 
AC167303.2) was obtained, which contains all introns and 
exons of the porcine IRF6 gene, and then the genomic 
organization was also subsequently deduced by 
experimental confirmation using exon-exon PCR across 
predicted introns or exon-intron PCR. Finally, the sequence 
of porcine IRF6 genomic DNA was submitted to NCBI 
GenBank database. 

The porcine IRF6 gene consist of seven exons, using 
porcine genomic DNA, seven primer pairs were designed 
for PCR amplification of DNA sequences including all exon 
regions of the IRF6 gene (Table 1). The DNA from two pigs 
from each breed was selected for PCR ampUfication. All 
PCR fragments were purified with a Gel Extraction Mini 
Kit (Beijing Tiangen Biotechnology, China) and then 
sequenced. After comparing the sequences of obtained from 
the same PCR fragments from two Landrace, two Large 
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Table 1. Primers used for mRNA expression and SNP detection of the porcine IRF6 gene 


Fragments 


Primers sequence (5 '-SO 


Product size (bp) 


PGR (Tm) 


Used for 


IRF6-RT 


CTGAAGCCCTGGCTGGTAG 
TGATGTTCAGGAAGGGGAAG 


436 


60.6 


mRNA expression 


IRF6-1 


TGAGGGAGCTGGAAAACAAC 

CAGATCATGGCCCTCCAC 


228 


60.2 


SNP detection 


IRF6-2 


TTCGTAAGCGAGGCAATTTT 
GGCAGTTTCCTGTCTTCCAG 


291 


59.6 


SNP identification 


IRF6-3 


TGATGTTCAGGAAGGGGAAG 
AGGCCTGTGCATCTGTTTCT 


152 


59.8 


SNP identification 


IRF6-4 


GCAGAAGGTGGACAGGAAAT 
GCTCACCCACATCCTTGTCT 


253 


59.5 


SNP identification 


IRF6-5 


CAGTTCCCACCATCCATCTT 
TTCAGTACCGTGGGAAGGAG 


410 


59.7 


SNP identification 


IRF6-6 


AAGGCCAAGAGATGTTGACG 

GATCTTATCGCCCACCAGAA 


156 


60.0 


SNP identification 


IRF6-7 


TGAAGTCATCAAGGCAGAGC 
GrGGCGCGAAlGArCTArG 


612 


59.5 


SNP detection 



White and two Songliao Black pigs using DNAMAN 
software (version 5.2.10, Lyn-non Biosoft, Canada), 
polymorphisms of IRF6 gene were found. Then all 302 
samples from three pig populations were genotyped by 
directly sequencing of PCR products respectively. 

Association analysis 

In the end, we performed an association test between the 
genotypes of SNPs and three immune traits by SAS 
software (version 9.13), based on the following mixed 
model: 

y = Xp+Zb+e 

Where y is the vector of phenotype for immune traits 
analyzed; P is the vector of fixed effects parameter 
including breed, ELISA plate effect and genotype effect, 
and X is its incidence matrix; b is the mixed vector of 
random component including sires effect and dam effect 
within sires, and Z is its incidence matrix; and e is the 
vector of residual effect. 



RESULTS AND DISSCUSSION 

Sequence analysis of porcine IRF6 gene 

The porcine IRF6 cDNA is 2,203 bp in length 
(GenBank Accession Number: AF327368) and contains a 
1,404 bp open reading frame (ORF) encoding a 467-amino 
acid polypeptide. Certainly, a polyadenylation stretch is at 
the 3' end of the sequence, while a putative polyadenylation 
signal AATAAA is at 1,765 to 1,771 bp from it. The 
complete porcine 1RF6 genomic DNA sequence spans 
13,720 bp (GenBank Accession Number: HQ026023) and it 
consists of 6 introns separated by 7 exons (Figure 1). The 
intron/exon organization of porcine (6 introns and 7 exons) 
differs from that of human IRF6 (8 introns and 9 exons) and 
mouse 1RF6 (7 introns and 8 exons) (Ben et al., 2005). 
These exons in coding region were also compared with that 
of human; but exon sizes of porcine IRF6 gene were not 
identical with that of human (Figure 1). All IRFs share a 
highly conserved helix-turn-heUx DNA-binding domain 
(DBD) and a less conserved protein binding domain 
(Taniguchi et al., 2001). The DBD is defined by 5 
tryptophan (W) residues which are separated by 10 to 18 



Ptg IRF6 



Humaa IRF6 5'- 



ExoBl 



98 



BkoB2 £xoii3 Enm4 ExonS Enmfi Exob? 



Exoal Ex(n2 ExonS Exoa4 ExonS Exob6 Exob7 ExoaS ExonP 



Md.»IRF6 .-• 1||3916 177 4330 205^84 129 . n.. 159 i,,4^...119...^76y ^ 



Exonl Exou2 Exon3 Exon4 ExonS Exon6 £xon7 ExodS 



Figure 1. Genomic organization comparisons among pig, human and mouse IRF6 gene. Exon sequences are represented in the boxed 
areas (length in bases on top). 
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amino acids respectively (Escalante et al., 1998; Liu et al., 
2011). Through the DBD, IRFs bind to a core DNA 
sequence GAAA, termed interferon stimulating response 
element (ISRE) and IFN regulatory element (IRF-E) 
(Tanaka et al., 1993). The porcine IRF6 displayed a higher 
degree of homology than other known IRFs and also 
contain 5 tryptophan residues that are each separated by 10- 
18 amino acids in the first 100 amino acids of the amino- 
terminal region, as weU as a high sequence homology in 
their N-terminal region in DBD. In contrast, they showed 
very low homology in the C-terminal region. These 
sequence characters could contribute to further 
understanding of its functions in viral defense and immune 
regulation. 

Tissues expression of porcine IRF6 gene 

To determine tissue distribution of the porcine IRF6 
gene, reverse transcription quantitative PCR (RT-PCR) was 
applied to nine different tissue samples. Results showed that 
the mRNA of IRF6 was expressed in all analyzed tissues 
except spleen (Figure 2). The expression levels of IRF6 
detected in lung, hver and lymph node were higher, 
comparatively, while in brain and heart the levels were 
lower. We did not find any expression of the IRF6 gene in 
the spleen. In some previous studies, 1RF6 gene of Xenopus 
was found to be expressed as a maternal transcript and later 
in the posterior somite mesoderm of neural groove stage 
embryos, and IRF6 expression has been observed only in 
the ovary after adulthood was attained (Hatada et al., 1997). 
Moreover, Kondo et al. (2002) detected the expression of 
the mouse IRF2 mRNA in several tissues by RT-PCR and 
also found higher expression levels in lung, liver and much 
lower expression levels in brain, heart, which is consistent 
with our results in pigs. 

Polymorphisms detection of porcine IRF6 gene 

Sequence comparisons of the same PCR fragment 
among three pig populations detected one SNP 
(HQ026023:1383 G>C) in exon? and two SNPs 
(HQ026023:130 G>A; 232 C>T) in the 3 promoter region 
of the porcine IRF6 gene. The three SNPs are all 
synonymous mutations, and do not induce a substitution of 



M Hb Ki Sp SM LN Lu Li Th Br 




Figure 2. mRNA expression of porcine 1RF6 gene was detected in 
different tissues by RT-PCR. He, heart; Ki, kidney; Sp, spleen; 
SM, skeletal muscle; LN, lymph node; Lu, lung; Li, Uver; Th, 
thymus; Br, brain; M, 100 bp DNA marker. 

amino acids. The three SNPs were genotyped by directly 
sequencing of the PCR products of all 302 pigs. The 
genotypes and allele frequencies of the identified SNPs in 
three different pig populations were presented in Table 2. 
For two SNPs in 5' promoter region (HQ026023:130 G>A 
and HQ026023:232 OT), the allele G and C have higher 
frequencies than allele A and T in Large White and Songliao 
Black pig populations, but have lower allele frequencies in 
Landrace population. In addition, for the SNP 
((HQ026023:1383 G>C), the allele G had higher allele 
frequencies, while there was only genotype GG in Chinese 
indigenous breed Songhao Black pig population. 

Association analysis of porcine IRF6 gene in three pig 
populations 

Based on the aforementioned work, we performed 
association analysis to confirm whether the identified SNPs 
in the breeds have a relationship with immune traits. 
Results showed that all the traits among three breeds were 
significantly different except for the ratio of IFN-y to ILIO 
(d 35) in serum (p<0.05) (Table 3). The serum IFN-y level 
(d 20; d 35) of Landrace was significantly higher than Large 
White and Songliao Black (Table 3); In contrast, Songliao 
Black had a significantly higher serum ILIO level (day 20; 
day 35) than Large White and Landrace (p<0.01) (Table 3); 
The ratio of IFN-y to ILIO (d 20) in serum of Songliao 
Black was significantly lower than Landrace and Large 
White (p<0.05) (Table 3). As reported in a previous study, 
the capacity to activate or inhibit monocytic and T 
lymphocytic functions could be reflected by the ratio of 



Table 2. Genotype and allele frequencies of the SNPs detected in the porcine 1RF6 gene' 

HQ026023:g.l30 G>A HQ026023:g. 232 C>T HQ026023:g. 1383 G>C 



Breed Genotype Allele Genotj'pe Allele Genotype Allele 

frequencies frequencies frequencies frequencies frequencies frequencies 





GG 


GA 


AA 


G 


A 


cc 


CT 


TT 


C 


T 


GG 


GC 


cc 


G 


C 


Landrace (68) 


0.25 


0.37 


0.38 


0.44 


0.56 


0.25 


0.37 


0.38 


0.44 


0.56 


0.84 


0.16 


0 


0.92 


0.08 


Large White (163) 


0.31 


0.44 


0.25 


0.53 


0.47 


0.28 


0.47 


0.25 


0.52 


0.48 


0.40 


0.41 


0.19 


0.61 


0.39 


Songliao Black (71) 


0.32 


0.58 


0.10 


0.61 


0.39 


0.31 


0.56 


0.13 


0.59 


0.41 


1.00 


0.00 


0.00 


1.00 


0.00 



The number of animals in each breed is indicated between parentheses. 
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Table 3. Association analysis of immune traits in three pig breeds 






Breeds (Least square mean+SE) 




Traits 


Landrace 


Large White 


Songliao Black 




(n = 68) 


(n = 163) 


(n = 71) 


IFN-Y (d 20) 


107.57±12.61^ 


58.29112.98" 


59.86113.85'' 


IFN-y (d 35) 


100.47+15.64"^ 


59.52+16.10'' 


66.50+17.17" 


IFN-y/ILlO (d 20) 


1.53+0.47' 


1.32+0.49" 


0.66+0.52" 


IL10(d 20) 


133.62±41.29^ 


81.82±42.49^ 


173.74145.32" 


IL10(d35) 


100.39±33.50'' 


54.38134.48" 


107.21136.78" 


IFN-y/ILlO (d 35) 


2.152+0.45 


1.91+0.46 


1.6510.50 



' '' Statistically different of least square means (p<0.05). '^^ Statistically different of least square means (p<0.01). 



IFN-y/IL-lO production, and a higher ratio has also been 
shown to have relevance with depressive disorders (Maes, 
1999). In our study, the ratio of IFN-y to IL-10 in Songliao 
Black pig was lower, which could be explained as the 
Chinese indigenous breed has a better T lymphocytic 
balance and immune capacity than other two western 
commercial pig breeds in our study. Besides, the results also 
indicate that genetic background is one of the important 
factors that influence the immune capacity. 

On the other hand, the results of association analysis 
between SNPs and immune traits also showed that the SNP 
(HQ026023:1383 G>C) of IRF6 gene was significantly 



associated with the level of IFN-y (d 20) in serum (p = 
0.038) and the ratio of IFN-y to ILIO (d 20) in serum (p = 
0.041) (Table 4). The other two SNPs (HQ026023:130 
G>A; 232 C>T) were highly significantly associated with 
ILIO level in serum both at the d 20 (p = 0.005; p = 0.001) 
and the day 35 (p = 0.004; p = 0.006) (Table 4). In the 
further analysis, for the SNP (HQ026023:1383 G>C), the 
IFN-y level in serum of pigs associated with genotype GG 
were significantly higher than with genotype CC (p<0.05); 
but the ratio of IFN-y to ILIO with genotype CC was 
significantly higher than those with genotype GG (p<0.05). 
For the other two SNP (HQ026023:130 G>A; 232 C>T), 



Table 4. Association analysis and multiple tests of the IRF6 gene SNPs with immune traits in three pig populations 



p-value 



Immune traits HQ026023 HQ026023 HQ026023 

130 G>A 232 OT 1383 G>C 



Polymorphism genotypes 
(Least square mean+SE) 



IFN-y 0.190 0.216 0.038* HQ026023:g.l383 G>C 

(d20) GO GC CC 

44.59+26.18' 34.08+26.45" 33.87+26.64" 

IFN-y 0.841 0.941 0.389 

(d 35) 

IFN-y/ILlO 0.616 0.703 0.041* HQ026023:g.l383 G>C 

(d20) GG GC CC 

0.89+0.56 " 0.9010.47" 1.08+0.34" 
ILIO 0.005** 0.001** 0.317 HQ026023:g.l30 G>A 

(d20) GG GA AA 

349.58+52.54^ 200.12170.06" 138.93+99.99" 
HQ026023:g.232C>T 

CC CT TT 

76.89+57.65"^ 81.44+69.56"^ 406.23+95.59" 
ILIO 0.004** 0.006** 0.819 HQ026023:g.l30 G>A 

(d35) GG GA AA 

279.01+42.50"^ 162.64+56.67" 72.33+80.88" 
HQ026023:g.232 OT 

CC CT TT 

12.29+46.63"^ 55.41+56.27"^ 326.19+77.32" 

IFN-y/ILlO 0.280 0.320 0.643 

(d 35) 



* p<0.05; ** p<0.01. " Statistically different of least square means (p<0.05). *^ Statistically different of least square means (p<0.01). 
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the ILIO level in serum of pigs associated with genotype 
GG and TT were significantly higher than those with 
genotype AA and CC (p<0.01). 

IRF6 plays an important role in initiating cellular 
interferon- stimulated gene-mediated antiviral responses and 
in the transforming growth factor-P signahng pathway. 
Moreover, in human, this gene is essential for the 
development of the lip and palate and involved in the 
development of the skin and external genitalia. Mutations in 
IRF6 gene are pathogenic for both Van der Woud 
syndromes (VWS) and Popliteal pterygium syndrome (PPS) 
(Kondo et al., 2002). Until now, two missense mutations in 
the IRF6 gene had been identified in two Japanese families 
with PPS (Matsuzawa et al., 2010). Another two novel 
mutations which may affect splicing in the IRF6 gene 
associated with VWS were also reported by Scioletti et al. 
(2010). In the pig, another IFN regulatory factor (IRFl) has 
been reported to be significantly associated with cytokine 
traits (Liu et al., 2010; 2011). The SNP (HQ026023:1383 
G>C) in exon7 and two SNPs (HQ026023:130 G>A; 232 
C>T) in 5' promoter region of IRF6 were not missense 
mutations or induced amino acid substitution. However, the 
synonymous SNPs can also affect protein expression by 
alteration or by increasing the stability of the mRNA 
(Capon et al., 2004). Additionally, a "silent" polymorphism 
changes substrate specificity (Kimchi-Sarfaty et al., 2007). 
Furthermore, it is known that mutations in the promoter 
region are involved in biological processes such as 
transcriptional and posttranscriptional control of gene 
expression (Nakaya et al., 2007). 

As the immune system plays an essential role in disease 
resistance of animals, genes involved in the response of 
immune system could be regarded as the candidate genes 
for selection. Our results provided a straightforward insight 
that the 1RF6 gene which has effects on serum IFN-y, lL-10 
level and ratio of IFN-y to IL-10 in serum could served as a 
genetic marker for inmiune traits. But the number of pigs 
analyzed in our study was restricted, further investigations 
are needed to confirm the relationships between the SNPs 
with immune traits among other pig populations. 

In summiiry, the 13,720 bp-length IRF6 genomic 
organization with seven exons and six introns was identified. 
The result of mRNA expression revealed that the porcine 
IRF6 gene was expressed widely in lung, liver, lymph node, 
brain, kidney, skeletal muscle, thymus and heart but not in 
the spleen. One SNP (HQ026023:1383 G>C) in exon 7 and 
two SNPs (HQ026023:130 G>A; 232 C>T) in 5' promoter 
region of IRF6 gene were found. Further association 
analysis showed that the three SNPs were highly 
significantly associated with serum IFN-y (d 20), ratio of 
IFN-y/ILlO (d 20) and ILIO level (d 20; d 35) respectively. 
The results indicated that IRF6 gene could be regarded as a 



genetic marker with effects on these immune traits in 
porcine disease resistance breeding. 
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